The turbulent mean-velocity profile: it is all in the spectrum 
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It has long been surmised that the mean-velocity profile (MVP) of a pipe flow is closely related 
to the spectrum of turbulent energy. Here we perform a spectral analysis to identify the eddies that 
dominate the production of shear stress via momentum transfer. This analysis allows us to express 
the MVP as a functional of the spectrum. Each part of the MVP relates to a specific spectral 
range: the buffer layer to the dissipative range, the log layer to the inertial range, and the wake 
to the energetic range. The parameters of the spectrum set the thickness of the viscous layer, the 
amplitude of the buffer layer, and the amplitude of the wake. 



Although most flows in nature and technology are tur- 
bulent flows over confining walls, these flows have re- 
mained amongst the least understood phenomena of clas- 
sical physics [1| . Take the simplest example (an example 
whose applications are legion): the turbulent flow in a 
long cylindrical pipe with a cross section of radius R and 
a smooth internal wall. If the flux is kept steady, the 
velocity of the flow at a distance y from the wall of the 
pipe may be averaged over a long period of time to ob- 
tain a mean velocity u. The function u{y) is called the 
mean- velocity proflle (MVP), and there is a MVP for 
each value of the Reynolds number. Re = UR/v, where 
U is the mean velocity of the flow (i.e., the flux divided 
by the cross-sectional area of the pipe) and v is the kine- 
matic viscosity of the fluid 0, Q • (Re quantifles the rel- 
ative importance of inertia and viscosity in the flow; the 
higher the Re, the more turbulent the flow.) MVPs were 
first measured 70 years ago [3| and have recently been the 
subject of exacting experiments and numerous compu- 
tational simulations [y, 0] ■ Theory has meanwhile lagged 
well behind experiments and simulations. 

The first theory came soon after the earliest exper- 
iments. Ludwig Prandtl showed that when plotted in 
terms of the dimensionless "wall variables" u and y, the 
MVPs for different values of Re collapse into a single 
MVP close to the wall (Fig . [J) 3] . Th e wal l variables 
are defined by u = u/ \Jt^I p and y = {^Jtq/ p)/v, where 
p is the density of the fiuid and ro is the shear stress — or 
shear force per unit area — that develops between the fiow 
and the wall of the pipe. Prandlt also showed that over 
most of its domain u{y) follows the "log law of the wall," 



u{y) = -\ny + B, 

K 



(1) 



where k (the "Karman constant" ) and B are dimension- 
less constants which can be estimated by fitting experi- 
mental data [1]. 

Numerous variants of the log law of the wall have been 
proposed since the time of Prandtl, and recently it has 
been argued that the MVPs for different values of Re 
might not be logarithmic after all. Instead, they might 
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FIG. 1: Log-linear plots of the MVPs in the "wall variables" u 
and y for four values of Re. The symbols are from experiments 
[Sj and the solid lines from simulations 0] • Each MVP extends 
from the wall (which corresponds always to y — 0) to the 
centerline of the pipe (which corresponds to a value of y that 
depends on Re); the MVPs collapse on a single MVP close 
to the wall. Focusing on any one MVP, we can parse the 
horizontal axis from left to right to find the "viscous layer" 
(where the MVP has a positive curvature), the "buffer layer" 
(where the MVP has a noticeable negative curvature), the 
"log layer" (where the MVP follows the log law of the wall 
((T])), and a "wake" (where the MVP overshoots the log law of 
the wall close to the centerline of the pipe) @| . 



be power laws with Re-dependent exponents and a loga- 
rithmic envelope . Yet the original theory of Prandtl 
and this recent theory of Barenblatt have been predicated 
on dimensional analysis and similarity assumptions (they 
differ in the similarity assumptions 9]), without reference 
to the turbulent spectrum. As a result, these theories 
cannot be used to relate, for example, the log law of the 
wall or the constants thereof to the turbulent spectrum. 
(The turbulent spectrum may be said to embody the fab- 
ric of a turbulent state; in practice it is an account of the 
form in which the kinetic energy is apportioned among 
the eddies of different sizes in the flow.) Our aim here 
is to find the missing link between the MVP and the 
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turbulent spectrum. 

We adopt the phenomenological imagery of "turbu- 
lent eddies" [10, Kill and use the spectrum of tur- 
bulent energy II4I at the wavenumber fc, E{k), 
to determine the velocity of the eddies of size s, 



in the form 



l{Kesf/^k^^/^Cd{rik)ce{Rk). Here 
less parameter, e is the turbulent power per unit mass 



where E{k) — 
is a dimension- 



,3/4^-1/4 



is the viscous lengthscale 2], 
2 



R is 



the largest lengthscale in the flow, | (nge)^^^ k ^^/^ is the 



Kolmogorov spectrum IJ, [15| , and Cd and Ce are dimen- 
sionless correction functions — the dissipative-range cor- 
rection and the energetic-range correction, respectively. 
For the dissipative-range correction we adopt the usual 
exponential form, Cd{r]k) — exp(— /S^f^A;), and for the 
energetic-range correction the form proposed by Karman, 

are non- 
negative dimensionless parameters 
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2 roo , 
s in the inertial range 



By introducing the dimensionless variable ^ = sfc, we 
can write Vg = {n^esY/'^y/l , where I = X{r]/s, s/R) = 
J^^-5/3exp{-^(3dv/s) (l+/3e(s/ii)Ve')-''/'rfe For 
<^ s <^ R), I — 1, and therefore 
with the implication that the 
velocity of an eddy of the inertial range (i.e., an eddy 
of size 77 ^ s <C i?) increases with the size of the eddy. 
The same result, Vs = {k^ss)^^^, may be derived from 
Kolmogorov's four-fifth law, which also gives the estimate 
Ke = 4/5 16]. For s outside of the inertial range, / < 1. 
It follows that an eddy of the dissipative range or the 
energetic range (i.e., an eddy of size s » rj or s ~ R, 
respectively) has a velocity < (K^es)-^/'^ — i.e., the eddy 
is slower than an imaginary eddy of the same size in the 
inertial range. 

Consider now the flow in a smooth-walled pipe. The 
energy equation gives and expression for e in the form 
e = Ttu'/p, where is the turbulent shear stress, u' = 
du/dy, and y is the distance to the wall [l^.The sum of 
the turbulent shear stress Tt and the viscous shear stress 
pvu' is the total shear stress r, which in a pipe of radius 
R is given by r = To{l — y/R), where tq is the shear stress 
at the wall [I4I . Using the definition of the friction factor, 
/ = To/pU'^, we can write tq = pfU"^ [III; it follows that 



Tt + piyu' = pfU''{l~y/R) 



and 



fU\l 



y/R)u' 



(2) 



(3) 



We now seek to derive an expression for the turbulent 
shear stress Tt. Let us call Wy the wetted surface at a dis- 
tance y from the wall of the pipe (Fig. [2]). The turbulent 
shear stress that acts on Wy is produced by eddies that 
straddle Wy and transfer momentum across Wy (Fig. [2]). 
Thus an eddy of size s carries fluid of high horizontal 
momentum per unit volume (about pu{y + s/2)) down- 
ward across Wy and fluid of low horizontal momentum 
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FIG. 2: Schematic for the derivation of the turbulent shear 
stress. u{y -\- s/2) is the mean velocity at a distance y -f s/2 
from the wall; u{y — s/2) is the mean velocity at a distance 
y — s/2 from the wall. 



per unit volume (about pu{y— s/2)) upwards across Wy, 
and the eddy may be said to span a momentum contrast 
p{u{y + s/2) - u{y + s/2)) « psu'{y) The rate of 

momentum transfer across Wy is set by the velocity of 
the eddy, Vg (i.e., the velocity normal to Wy). Therefore, 
the turbulent shear stress produced by an eddy of size 
s scales as the product of the momentum contrast times 
the rate of momentum transfer, or psu'{y) Vg. 

In order to identify the dominant eddies that strad- 
dle Wy (i.e., the eddies that dominate the momentum 
transfer across Wy), recall that Vs (and therefore svg) in- 
creases with s. It follows that the larger the eddy the 
larger the turbulent shear stress produced by the eddy. 
Nonetheless, eddies much larger than y do not properly 
straddle Wy and can provide only a negligible velocity 
normal to Wy, and therefore a negligible turbulent shear 
stress. (This observation is purely a matter of geome- 
try.) We conclude that the dominant eddies that strad- 
dle Wy are the eddies of size s = y. The turbulent shear 
stress at a distance y from the wall is thus given by 
Tt = Krpyvyu'{y), where Kf IS cL dimensionless propor- 
tionality constant. 

We are now ready to obtain an equation for the MVP. 
We substitute ^ invy — {Keey)^^'^\/l, substitute the re- 
sulting expression for Vy in Tt = KT-pyVyu', and substitute 
the resulting expression for Tt in ((2]). After some algebra, 
we write the final result in terms of Re ^UR/v and the 
dimensionless variables y = y/R and u = u/U: 



K^rfu'^ + Rc^^u' - /(I -y) = 0. 



(4) 



where k = [k^k^)^/'^, u' = du/dy and / — I{rj/y,y). 
To obtain an equation for rj/y we substitute ([3]) in ry = 
j^3/4g-i/4 g^^^ change variables to y and ii: 



rj/y ^ Re-1/2 (/Reu'(l _ y) - u'2)-l/4 fj-\ 



(5) 



If for a fixed Re we let y — > 0, then rj/y ^ 00 (from 
([5])), / l{oo,y) = 0, and (|4]) simplifies to u' — /Re, 
which is the law of the viscous layer. If for a fixed 
y ^ 1 we let Re 00, then rj/y (from ^), y is 
in the inertial range (where / = 1), and ^ simplifies to 
it' — y/J / Ky, which we recognize as the log-law of the 
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wall with a Karman constant k = {k^k^ Note that 
the log law of the wall prevails where y is in the iner- 
tial range; it follows that the dominant eddies in the log 
layer are eddies of the inertial range. The presence of 
and Kr in the expression, k = (k^k^)^/^, reminds us of 
the underpinnings of the theory: the spectrum and the 
momentum transfer, respectively. As is fixed by Kol- 
mogorov's four-fifth law, which gives = 4/5 [la. Iist. 
the Karman constant k is but an alternative form of the 
momentum-transfer constant Kt- From k= 0.42 (the ex- 
perimental value 0), we have Kr = 0.34. 

The law of the viscous layer and the the log law of 
the wall may be made invariant to changes in Re and 
/, in the form u' = 1 and u' — l/ny, by choosing y = 
^^V7y = ^'i\fly/R and u = u/\fj — u/U\/J, which 
we recognize as the wall variables. In the wall variables 
(|4]) becomes 

^2/3^2-/2 + ^' _ (1 _ ^/Re V7) = 0, (6) 
where / = T{r]/y,y/Ke\/f), and ([5]) becomes 

7j/y = {u\l - y/Re V7) - u'^r'^^ F (7) 

From these equations it is apparent that in the wall vari- 
ables there is a single MVP except close to the centerline 
of the pipe, where y « Re-\/7 (or y R). 

We now ascertain under what conditions (O and ([7]) 
are compatible with a nonvanishing turbulent shear stress 
close to the wall. Suppose that at a point y <C ReyT (or 
?/ <C i?) we have / = X{ri/y, 0) > 0, and therefore n > 0. 
Then, we can eliminate u' from ([6]) and ([7|) to obtain 

, ^ f {v/y)^/^ + ^^/'lHv/y,0) \ ^'^ 

V \ nV^^/y)mi{j^/y^Q) ) ■ ^ ^ 

Plots of y vs. rj/y (Fig. [3]) reveal that for any given k and 
Pd there exists a minimum value of y, to be denoted y^. 
Therefore, for y < jjy it must be that / — 0, rt = 0, and 
u' — 1. As this latter equation is the law of the viscous 
layer, we identify jjy with the thickness of the viscous 
layer [§]. Note that ijy depends on the dissipative-range 
parameter (3d (Fig. [Sb), and that for Pd = there is no 
viscous layer (?/„ = 0). 

We are now ready to compute the MVP. For simplicity 
we use the Blasius relation for the friction factor, / = 
0.033 Re~^^^ 3. Values of the spectral parameters Pd 
and f3e from the fitting of spectra are of 0(1) [12]; here 
we use /3d = 1 and (3e — 4. The thickness of the viscous 
layer is set by f3d and k; for /3d = 1 and k = 0.42, ?/.„ = 3.9 
(Fig. [3]d). Therefore, for y < 3.9 we can write uijj) = y, 
and for y > 3.9 we must compute u{y) by integrating ([6]) 
with boundary condition {((3.9) — 3.9. The results are 
shown in Fig. Given the spectrum, the theory yields 
the entire MVP with all of its distinctive features. The 
specific conection between each one of these features and 
the spectrum will become apparent in what follows. 
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FIG. 3: (a) Plots of (O for «: = 0.42 and three values of P^. 
(b) Plots of the thickness of the viscous layer y„ as a function 
of Pd for three values of k. 



To elucidate the effect of the energetic-range correc- 
tion, we recompute the MVP using /3e = 7 (a larger value 
than before) and Pe — ^ (the smallest possible value, 
which corresponds to having no energetic-range correc- 
tion). The results (Fig. (Hd) indicate that the energetic- 
range correction steepens the MVP in the wake. The 
dominant eddies in the wake must therefore be eddies of 
the energetic range. These eddies are slowed down by the 
energetic-range correction (and made less adept at trans- 
fering momentum). This effect explains the steepening 
of the MVP in the wake. 

We have seen that the dissipative-range correction sets 
the thickness of the viscous layer (Fig. [SId)). To under- 
stand further the effects of the dissipative-range correc- 
tion, we recompute the MVP using a few values of Pd, 
including Pd = Q (the smallest possible value, which cor- 
responds to having no dissipative-range correction) . The 
results (Fig. |4):) indicate that the dissipative-range cor- 
rection causes the buffer layer to form, so that for Pd = ^ 
there is no buffer layer. (The buffer layer is the part of 
the MVP where the MVP has a negative curvature |8j]-) 
The dominant eddies in the buffer layer must therefore be 
eddies of the dissipative range. These eddies are slowed 
down by the dissipative-range correction; the larger the 
eddies, the less they are slowed down, and the more adept 
they remain at transferring momentum. As the size of the 
dominant eddies increases with the distance to the wall, 
the MVP becomes less steep as we traverse the buffer 
layer from the outer edge of the viscous layer (where the 
eddies are fully viscous) to the inner edge of the log layer 
(where the eddies are fully inertial). This effect explains 
the negative curvature of the MVP in the buffer layer. 

It is apparent from Fig. Uj; that the constant B of the 
log law of the wall is set hy Pd- A plot of S as a function 
of Pd is shown in Fig. [4ji. 

Interestingly, the MVP for Pd — (^ (no dissipative-range 
correction) can be obtained analytically everywhere save 
the wake. In fact, for /3d = and y <C Rev7 (or y <^ R), 
I — 1 and ^ simplifies to K^y^^'S _j_ _ _ which 
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FIG. 4: (a) The MVPs computed using k = 0.42, /3e = 3, and 
l3d ~ 1- The thick dot indicates the point of contact between 
the viscous layer and the buffer layer, (b) The same as in 
(a) but using (3^—0 (dashed lines) and /3e = 7 (solid lines). 
Inset: detail of the wakes, (c) The same as in (a) but using 
Pd = (bottom), 1, 2, and 3 (top). Inset: A plot of Q for 
K = 0.42. (d) A plot of B as a function f3d for k = 0.42. 



with u(0) = yields 



arcsinh(2Ky) 1 — ^1 + K^y^ 



(9) 



This is what the MVP would be away from the wake 
if the Kolmogorov spectrum were valid even for vanish- 
ingly small eddies (inset of Fig. For <C <C Rey^ 
(or ^ y ^ i?), (m simplifies asymptotically to w '--^ 
(1/K)lny + B, withS = (-l + ln4K)/K. Thus, for /3d = 
and K = 0.42, B = —1.15 in accord with Fig.HH. 

We have established the long-surmised link [3| between 
the mean velocity profile and the turbulent spectrum. To 
test our results, we have shown that the usual model of 
the spectrum (a power-law inertial range with correc- 
tions for the dissipative range and the energetic range) 
is in itself sufficient to compute with no additional as- 
sumptions a mean velocity profile complete with viscous 
layer, buffer layer, log layer, and wake. The thickness 
of the viscous layer, the two constants of the log law of 
the wall, and the amplitude of the wake are all set by 
the dimensionless momentum-transfer constant and 
the usual spectral parameters — the parameter /3e of the 
energetic-range correction and the parameter f3d of the 
dissipative-range correction. The relation between a spe- 
cific feature of the MVP and the spectral parameters re- 
minds us of the underlying physics. Thus, for example, 
the Karman constant is independent of both (3e and (3d, 
and is therefore unaffected by the energetic-range and 
dissipative-range corrections, with the implication that 
the eddies that dominate the momentum transfer in the 



log layer are eddies of the inertial range. More broadly, 
the close relation between the mean velocity profile and 
the turbulent spectrum indicates that in turbulence, as 
in continuous phase transitions, global variables are gov- 
erned by the statistics of the fluctuations. 
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